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Abstract. We use quantum Monte Carlo methods in the framework of the interacting nuclear shell model to calculate the 
statistical properties of nuclei at finite temperature and/or excitation energies. With this approach we can carry out realistic 
calculations in much larger configuration spaces than are possible by conventional methods. A major application of the 
methods has been the microscopic calculation of nuclear partition functions and level densities, taking into account both 
correlations and shell effects. Our results for nuclei in the mass region A ~ 50 — 70 are in remarkably good agreement with 
experimental level densities without any adjustable parameters and are an improvement over empirical formulas. We have 
recently extended the shell model theory of level statistics to higher temperatures, including continuum effects. We have also 
constructed simple statistical models to explain the dependence of the microscopically calculated level densities on good 
quantum numbers such as parity. Thermal signatures of pairing correlations are identified through odd-even effects in the heat 
capacity. 



1. INTRODUCTION 

The statistical properties of nuclei at finite temperature or 
excitation energy are important for nuclear astrophysics. 
Level densities are needed for estimates of reaction rates 
in nucleosynthesis jlll- Such reactions include in partic- 
ular neutron capture in the s and r processes and proton 
capture in the rp process. Nuclear partition functions are 
used in the calculation of thermal stellar thermal reaction 
rates [^. 

Also of interest are the signatures of phase transitions 
(e.g., pairing) in finite systems. Strictly speaking, phase 
transitions occur only in the thermodynamic limit (of 
bulk systems). In finite systems, fluctuations around the 
mean field are important and they smooth the singulari- 
ties associated with the phase transitions. An interesting 
question is whether signatures of these transitions remain 
despite the fluctuations. 

The fundamental quantity for calculating statistical 
properties at temperature T is the partition function 



Z{T) = Tre-"/^ 



(1) 



where H is the nuclear many-body Hamiltonian. Such a 
Hamiltonian is described by the interacting shell model. 
The shell model has been successful for describing the 
properties of, e.g., iii-shell and low /;?-shell nuclei. In 
such nuclei conventional diagonalization in a complete 
major shell is possible. However, in medium-mass and 
heavy nuclei the dimensionality of the model space is 
often too large to allow for exact diagonalization. 



Mean-field approximations are tractable but are not 
always sufficient. Correlations beyond the mean field 
can be taken into account by considering fluctuations 
around the mean field solution. Formally, this can be ex- 
pressed by a mathematical transformation known as the 
Hubbard-Stratonovich (HS) transformation |3]- In gen- 
eral, the HS transformation requires an integration over 
a large number of fluctuating fields. Quantum Monte 
Carlo methods have been introduced to integrate over 
these fields. In the context of the nuclear shell model 
this approach is known as the shell model Monte Carlo 
(SMMC) method K,"?,"]. We briefly review the SMMC 
method (Section |5} and then discuss its applications for 
the calculations of statistical nuclear properties (Sections 
EUandEll. 



2. THE MONTE CARLO APPROACH 
2.1. Hubbard-Stratonovich transformation 

The SMMC method is based on the HS transformation 
Pij, which describes the Gibbs ensemble e^l^^ at inverse 
temperature j3 (j3 = 1/T) as a coherent superposition of 
one-body propagators U^ 



-flH 



[cy]GaUa 



(2) 



Ga is a Gaussian weight and each Ua describes the 
imaginary-time propagator of non- interacting nucleons 
moving in external time-dependent auxihary fields (7(t). 



The thermal expectation value of an observable O is 
given by 



{O) 



Tr((9e-^^) _ /^[(j]Ga(0)aTrf/„ 



Tr(e-/5«) jSi[a]GaTxUo 



(3) 



where {0)a = Tr{OUa)/^^Ua is the expectation value 
of the observable O evaluated for a sample a of the 
auxihary fields. 

The calculation of the integrands in Eq. Q reduces to 
matrix algebra in the single-particle space. The one-body 
propagator Ua is represented by an Ns x Ns matrix Ucy 
in the single-particle space {Ns is the number of single- 
particle orbitals). The grand-canonical trace of Ua in 
the many-particle fermionic space can then be calculated 
from 

Trf/<, = det(l+U„). (4) 

Eq. describes the grand-canonical partition function 
of non-interacting fermions in time-dependent external 
fields C7(t). Similarly, {0)a can be expressed in terms of 
the matrix Ufj using Wick's theorem. 

In finite nuclei it is important to consider a definite 
number of protons and neutrons. Therefore the traces in 
Eq. (|3} should be evaluated in the canonical ensemble. 
Canonical traces can be calculated by exact particle- 
number projection. 

At sufficiently high temperatures, only static configu- 
rations of the fields a are important, leading to the static 
path approximation (SPA) |7]. Of particular importance 
are large-amplitude fluctuations of the relevant order pa- 
rameters. For example, in the Landau theory of shape 
transitions, static fluctuations of the quadrupole deforma- 
tions have explained the observed temperature and spin 
dependence of the giant dipole resonance 1 8]. 

At low temperatures, it is necessary to take into ac- 
count all fluctuations including the quantal fluctuations 
which are described by time-dependent configurations 
of the fields a. This requires an integration over a very 
large number of variables and in practice can be done by 
Monte Carlo methods. 



2.2. Monte Carlo methods 

The multi-dimensional integral is evaluated exactly by 
Monte Carlo methods, in which the a fields are sampled 
according to the distribution 



Wa = Ga\i:rUa\. 

The observables are then estimated from 

Za{0)a^a 



(O) 



La^o 



(5) 



(6) 



where <I><j = Ti Ua/ \Ti Ua\ is the sign of the one-body 
partition function TiUa 

Often the sign ^Pa is not positive for some of the 
samples a. When the statistical uncertainty of the sign 
is larger than its average value, the method fails. This 
known as the Monte Carlo sign problem. Most effective 
nuclear interactions suffer from this sign problem at low 
temperatures. A practical solution to this sign problem is 
discussed in Ref. |5]. Furthermore, "good-sign" interac- 
tions can be constructed for realistic calculations of col- 
lective properties. Such interactions are used in the ap- 
plications described below. 

In general, the computational properties of the Monte 
Carlo approach scale more favorably as a function of the 
number of single-particle orbitals Ns, enabling calcula- 
tions in much larger configuration spaces than are possi- 
ble with conventional methods (the Monte Carlo method 
scales as A^^ compared with an exponential scaUng of a 
direct diagonalization). 



3. NUCLEAR LEVEL DENSITIES AND 
PARTITION FUNCTIONS 

The Fermi gas model ignores important correlations. In 
practice, correlations are taken into account through em- 
pirical modifications of the Fermi gas model. In partic- 
ular, good fits to the data are obtained using the back- 
shifted Bethe formula (BBF) ^ 

p(£,)-^a-'/4(£.-A)-VVV^^(^^, (7) 

where a is the single-particle level density parameter and 
A is the backshift parameter However, a and A have to 
be adjusted for each nucleus, and it is therefore difficult 
to predict the level density to an accuracy better than an 
order of magnitude. 

The interacting shell model includes both shell effects 
and residual interactions and is therefore a good micro- 
scopic model for the calculations of level densities. Trun- 
cations suitable for the description of low-lying states 
cannot be used at finite temperature and thus complete 
major shells must be included. This requires calculations 
in large model spaces and SMMC is a suitable approach. 
In Ref. 1 10] we have introduced a method to calculate 
level densities in the SMMC method. 



3.1. Thermodynamic approach 

The level density can be obtained as the inverse 
Laplace transform of the partition function. The average 
level density is found when this transform is evaluated in 
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FIGURE 1. The SMMC level densities of three iron isotopes 
(symbols). Also shown are the experimental level densities of 
^^Fe and ^^Fe (solid lines). 



the saddle point approximation 



P{E): 



1 



AE) 



^/IkT^C 



(8) 



where S{E) is the canonical entropy and C is the canoni- 
cal heat capacity. 

In the Monte Carlo approach, we calculate the thermal 
energy from E{j3) = (//) and integrate the thermody- 
namic relation —dlnZ/dji = E{I5) to find the partition 
function Z(j3). The entropy and heat capacity in (|8} are 
calculated from S{E) = lnZ + j3E and C = -j3^dE/dj3, 
respectively. 



3.2. Level densities in the iron region 

We have applied the SMMC approach to calculate 
level densities in the mass region A ^^ 50 — 70 ITOt fllll . 
The shell model space includes the complete fpgg/j 
shell. We have constructed a good-sign interaction that 
includes attractive monopole pairing and multipole- 
multipole interactions. The multipole interaction terms 
are based on a surface-peaked interaction whose strength 
is determined self-consistently. The quadrupole, oc- 
tupole and hexadecupole terms are then renormalized by 
factors of 2, 1.5 and 1, respectively. The pairing strength 
is determined from odd-even mass differences. 

We have found remarkably good agreement with ex- 
perimental level densities without any adjustable param- 
eters. An example is given in Fig. [2 which shows the 
level densities of several iron isotopes versus excitation 
energy. The SMMC level densities (symbols) are com- 
pared with the experimental level densities (solid lines). 
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FIGURE 2. The single-particle level density parameter a 
(left panels) and the backshift parameter A (right panels) versus 
mass number A for different isotopes. The SMMC results fsolid 
squares) are compared with experimental values (x's) 1131 and 
empirical formulas (solid lines) Ll3i1 . FromRef. Illll 



The microscopically calculated level densities are well 
described by the BBF 0. By fitting the SMMC level 
densities to the BBF, we can extract the parameters a and 
A. In Fig. 12 we show the dependence of these extracted 
parameters (solid squares) on mass number A for man- 
ganese, iron and cobalt isotopes and compare them with 
experimental values (x's) |J_2]. Our microscopic results 
are an improvement over empirical formulas (solid lines) 
llSll . The parameter a varies smoothly as a function of 
mass while A exhibits odd-even staggering because of 
pairing effects. 



3.3. Extending the theory to higher 
temperatures 

The SMMC calculations of Section 13.21 are realistic 
for temperatures below ^ 1.5 — 2 Me V. At higher tem- 
peratures one must treat larger model spaces. Extend- 
ing SMMC to larger spaces is however computationally 
time-consuming, and we have recently developed an ap- 
proximate method to extend the theory to higher temper- 
atures while using SMMC only in the truncated space 
114J. 

In the independent-particle model, it is possible to 
treat the complete space including all single-particle 
bound states as well as the continuum. In particular, the 
many -particle grand canonical partition function Z^.^ can 
be expressed in terms of the single-particle energies e„ij 



and scattering phase shifts 5/y(e) 



InZ 



•GC 



'sp 



Y,{2j+l)\l^Hl+e-P^''-'r^^] (9) 



de5p{e)\n[l+e-P 



(e-M)l 



where 5p (e) = tt^ ' Y.ij {2j + 1 )d5ij/de is the continuum 
contribution to the single-particle level density. In the 
calculations we have used a Woods-Saxon potential with 
the parametrization of Ref. llSll 

The canonical partition function Z^ at fixed particle- 
number A^ can be obtained in the saddle-point approxi- 
mation InZiv ~ InZ^*^ -Pl^N- \ In (2;r((AA?)2)), where 
{{AN)^) is the variance of the particle number fluctua- 
tion. 

In the presence of interactions, we combine the fully 
correlated SMMC partition in the truncated space with 
the independent-particle model partition in the full space 
through 



InZv = lnZ,,„- + lnZ,p - InZ,,,., 



(10) 



Here Z,, and Z^jr are the partition functions in the pres- 
ence of interactions in the full and truncated spaces, re- 
spectively. The subtraction of the last term on the r.h.s. 
of Eq. ( I10> avoids a double counting of the truncated de- 
grees of freedom. The partition function ( I10> includes 
correlations and should be realistic up to T ^ 4 MeV (at 
higher temperatures it is necessary to take into account 
the temperature dependence of the mean field). 

The logarithm of the excitation partition function Z' = 
ZeP^o (/t'q is the ground-state energy) is shown in Fig.|3] 
as a function of temperature. The extended partition 
Zy (solid squares) is compared with the SMMC parti- 
tion in the truncated space Z,,.,,- (open squares) and the 
independent-particle model partition Z^p (dashed line). 
The solid line is a fit of the extended partition function to 
the partition function associated with the BBF lll4ll . 

The corresponding extended level density for ^^Fe is 
shown in Fig. |4] (solid squares) and is compared with 
the level density in the truncated space (open squares). 
The extended level density is well described by the BBF 
(solid line) with fixed a and A up to T ~ 4 MeV. 



4. PROJECTED LEVEL DENSITIES 

It is often necessary to know the dependence of the level 
density and partition function on the good quantum num- 
bers such as parity, spin and isospin. In SMMC this can 
be achieved by introducing the appropriate projections in 
the HS transformation. 




FIGURE 3. The partition function Z' as a function of tem- 
perature for ^^Fe. The SMMC partition in the truncated space 
(open squares) is combined with the independent particle 
model partition (dashed line) to give the extended partition 
function (solid squares). The solid line is a fit to a partition 
function associated with the BBF. From Ref. ^41 . 
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FIGURE 4. The extended level density of ^''Fe (solid 
squares) as a function of excitation energy is compared with 
the SMMC level density in the truncated fpggi2 shell (open 
squares). The solid line is a fit of the extended level density to 
the BBF. From Ref. \JM 



4.1. Parity distribution 

We have calculated (in SMMC) the projected energies 
E± for even- or odd-parity states as a function of /3 and 
then applied the method of Section lTTI to find the even- 
and odd-parity level densities p± \ 10, 16]. 

Contrary to the assumption often used in nucleosyn- 
thesis calculations, we find that in some nuclei p+ 7^ p_ 
even at the neutron resonance energy. In Fig. [S] we show 
the SMMC ratio P-/p+ (symbols) versus excitation en- 
ergy for three nuclei in the mass region A ~ 50 — 70. The 
crossover from one dominating parity at low excitations 
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FIGURE 5. The ratio p^j p^ of odd- to even-parity level 
densities versus excitation energy £v for ^''Fe, Ni and Zn. 
The microscopic SMMC results are compared with the results 
of a simple statistical model (solid lines). From Ref. Ilal 



to equally likely parities at higher energies depends on 
the particular nucleus. We have introduced a simple sta- 
tistical model 1 16] to estimate the odd-to-even parity ra- 
tio. The model assumes that the quasi-particle states with 
parity n{n being the parity with the smaller occupation) 
are randomly populated. We find for the ratio of odd- to 
even-parity partition functions Z_(j3)/Z+(j3) = tanh/, 
where / is the occupation of quasi-particle states with 
parity n. To mimic effects of the quadrupole-quadrupole 
interaction, we use deformed quasi-particle states. 

The results of our model (solid lines in Fig.|3 are in 
good agreement with the SMMC results (symbols). 



4.2. Isospin distribution 

Isospin is approximately a good quantum number in 
nuclei. The isospin dependence of level densities can be 
calculated in SMMC by exact isospin projection \VT\. 
Isospin projection also allows us to take into account 
the proper isospin dependence of the nuclear interaction. 
This isospin dependence of the nuclear interaction can 
lead to significant corrections in the total level density of 
N ^Z nuclei. 



5. HEAT CAPACITY AND THE PAIRING 
TRANSITION 

The pairing interaction leads to superconductivity in bulk 
metals below a critical temperature Tc- This phase tran- 
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FIGURE 6. Right panels: the heat capacities of ^"Fe (top) 
and ^^Fe (bottom). The extended heat capacity (solid squares) 
and the SMMC heat capacity in the /pgg /2 shell (open squares) 
are compared with the heat capacity in the independent-particle 
model (dashed line). For ^''Fe we also the BCS heat capacity 
(dotted-dashed). Left panels: the experimental heat capacities 
of Yb and '^' Yb (solid lines) in comparison with the Fermi 
gas results (dashed lines) [20]. Notice the similarity in the shape 
of the heat capacity between the theoretical and experimental 
results. 



sition is described by the BCS theory and predicts a dis- 
continuity in the heat capacity at Tc ] 18]. The BCS the- 
ory is valid in the limit when the single-particle mean 
level spacing d is much smaller than the pairing gap A. 
However, in finite nuclei, d (within a shell) is comparable 
to A, and fluctuations around the mean-field solution are 
important. An interesting question is whether signatures 
of the pairing transition can be still be observed despite 
the large fluctuations. 

The heat capacity in SMMC is calculated as a nu- 
merical derivative of the thermal energy. Such calcu- 
lations lead to large statistical errors at low tempera- 
tures. We have introduced a novel method that takes 
into account correlated errors and reduces the statisti- 
cal error by almost an order of magnitude ]19]. Re- 
sults for ^"Fe and ^^Fe are shown in the right pan- 
els Fig. |6l The open squares are the truncated calcula- 
tions (in the fpg<)i2 model space) ] 19], while the solid 
squares describe the extended heat capacity ]14]. The 
heat capacity is suppressed in comparison with the BCS 
heat capacity (dotted-dashed line). However, in the even- 
even nucleus ^"Fe we still observe a 'bump' around 
T ~ 0.8 MeV when compared with the heat capacity of 



the independent-particle model (dashed line). This bump 
disappears almost entirely in the odd-even nucleus ^^Fe 
and is a clear signature of the pairing transition. 

The heat capacity was recently measured in rare-earth 
nuclei |20]. Results (solid lines) are shown in the left 
panels of Fig. |6j for the even-even nucleus ^^^Yb and 
the odd-even nucleus '^'Yb. These experimental heat 
capacities have similar shapes to our calculated heat 
capacities. 



6. CONCLUSIONS 

We have calculated statistical nuclear properties using 
the SMMC method. Statistical properties of particular 
importance in nuclear astrophysics are nuclear level den- 
sities and partition functions. In the Monte Carlo ap- 
proach, fully correlated calculations are possible within 
complete major shells. We have extended our method to 
higher temperature by combining the correlated calcu- 
lations in the truncated space with independent-particle 
model calculations in the full space, including the con- 
tinuum. 

Projected level densities at fixed values of the good 
quantum numbers can be calculated by incorporating 
exact projection techniques in the SSMC method. We 
have used such projection methods to calculate the parity 
and isospin distributions. 

In the finite nucleus, fluctuations beyond the mean 
field are important and smooth the singularities of the 
phase transitions. A particularly interesting example is 
the pairing transition. BCS theory predicts a discontinu- 
ity in the heat capacity at the transition temperature, but 
in the finite nucleus this signature is suppressed. Nev- 
ertheless, we do find a 'bump' in the heat capacity of 
even-even nuclei around the pairing transition temper- 
ature, while such a signature is not seen in odd-even 
nuclei. We therefore conclude that in the fluctuation- 
dominated regime, thermal pairing correlations are man- 
ifested through strong odd-even effects. 
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